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Abstract 



J(i7o) = min{ J(f2), il convex, il G Sad}, 



We focus here on the analysis of the regularity or singularity of solutions Oq to shape optimization problems 
among convex planar sets, namely: 

q 

(~| . where Sad is a set of 2-dimensional admissible shapes and J : Sad M is a shape functional. 

Our main goal is to obtain qualitative properties of these optimal shapes by using first and second order 
optimality conditions, including the infinite dimensional Lagrange multiplier due to the convexity constraint. We 
prove two types of results: 

i) under a suitable convexity property of the functional J, we prove that Hq is a M^^'^-set, p G [1, oo]. This 
<^ , result applies, for instance, with p ^ oo when the shape functional can be written as J{il) = + 

■ P{^), where R{il) = F{\^1\, Ef{^l), Ai(fi)) involves the area \n\, the Dirichlet energy or the first 

lO i eigenvalue of the Laplace-Dirichlet operator Ai (il), and P{^) is the perimeter of 57, 

. ii) under a suitable concavity assumption on the functional J, we prove that is a polygon. This result 

applies, for instance, when the functional is now written as J(0) = R{il) — P{^), with the same notations 
as above. 



Keywords: Shape optimization, convexity constraint, optimality conditions, regularity of free boundary. 

1 Introduction 

The goal of this paper is to develop general and systematic tools to prove the regularity or the singularity of optimal 
shapes in shape optimization problems among convex planar sets, namely problems like: 

min{J(r2), convex, Q.^Sad\, (1) 

where Sad is a set of admissible shapes among subsets of and J : Sad — )■ M is a shape functional. Our 
main objective is to obtain qualitative properties of optimal shapes by exploiting first and second order optimality 
conditions on ([T]) where the convexity constraint is included through appropriate infinite dimensional Lagrange 
multipliers. 
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Our approach is analytic in the sense that convex sets are represented through adequate pai^ametrizations and 
we work with the corresponding "shape functionals" defined on spaces of functions. In particular, we will use the 
classical polar coordinates representation of convex sets as follows: 

O„:=|(r,0)G [0,oo)xM; r<-iy|, (2) 

where u is a positive and 27r-periodic function, often called "gauge function of It is well-known that 

Qu is convex <;=^ u" + u > 0. (3) 
Thus, Problem (dJ may be transformed into the following: 

min {j{u) := J(J1„), u" + u>0, uG I'ad} , (4) 
where Tad is a space of 27r -periodic functions which will be chosen appropriately to represent Sad in O- 



We obtain two families of results depending on whether j is "of convex type" or "of concave type". In the first 
case, we prove regularity of the optimal shapes. In the second case, we prove that optimal shapes ai^e polygons. 

i) "Optimal shapes are regular": under a suitable convexity property on the "main part" of the functional j, we 
prove that any solution uq of dU) is VF^'^, which means that the curvature of dQo = diluo is an function 
whereas it is a priori only a measure: see Theorems I2.4[|2.6l and Corollary 12.71 To that end, we simply use the 
first optimality condition for the problem ([T]). 

The functionals under consideration here are of the form J(il) = + C(il), where r{u) := R{0,u) has 

an L^-derivative and C is like ^ below and satisfies a convexity condition. As a main example, we consider 
R{Q) = F{\Q,\, Ef{^}), Ai(r2)) which depends on the area \Q\, on the Dirichlet energy Ef{Q,) and/or on the 
first eigenvalue Ai(r2) of the Laplace operator on Q, (with Dirichlet boundary conditions), and C{Q.) = P{^) 
is its perimeter: see Section [l!2l In this case, we actually prove that the optimal shape is W'^'°° which means 
that the curvature is bounded. 

ii) "Optimal shapes are polygons": next, we prove that, under a suitable concavity assumption on the functional 
j, for any solution uq of dUl, uq + Uq is (locally) a finite sum of Dirac masses, so that is (locally) 
a polygon: see Theorems 12.91 12.121 and Corollary 12.131 The proof of this result is based on the second 
order optimality condition for the problem ([T]). We apply this result to shape optimization problems where 
J(f]) = R{Q) — P{^) where = F{\Q.\, Ef{Q.), Ai(f])) with the same notations as above, see Section 
14.21 This application involves some sharp estimates on the second shape derivative of the energy which are 
interesting for themselves: see Section l4.3.2l 

Our examples enlighten and exploit the fact that, in the context of shape optimization under convexity constraint, 
the perimeter is "stronger" than usual energies involving PDE, in terms of the influence on the qualitative properties 
of optimal shapes: if it appears in the energy as a positive term, it has a smoothing effect on optimal shapes, and on 
the opposite as a negative term, it leads to polygonal optimal shapes. 

Dual parametrization: Since our results are stated for the analytic functionals dH), we may apply them to the dual 
parametrization of convex sets instead of the parametrization with the gauge function: each convex shape can also 
be associated to its support function hQ{9) = max{x • e*^, x G Q,},9 G T and ([T]) again leads to the problem: 

min h" + h>0, he Kd} , (5) 
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where j{h) := J{n^), f]^ being now the set whose support function is h, and Tad are all support functions of 
admissible shapes G Sad- In this framework, if j satisfies the suitable convexity property, the regularity result (i) 
above holds for Hq minimizer of ([5]). However, this regularity does not imply that the corresponding optimal shape 
Qq := ri'^" is regular, but it exactly means that Qq is strictly convex: see Section [34l 

The situation is more similar to the gauge representation when exploiting the results (ii). Indeed, when they apply, 
they imply that the optimal shape is polygonal as well: see Remai'k l2.14[ 

Situation with respect to previous results: The second family of results (ii) is an extension of previous results 
obtained in |[T4l by the two first authors for the specific following functional of "local type": 



where G = G{9, u,q) : T x [0, oo) x R — )• M is strictly concave in q. Among these functional, we find for 
instance the area the perimeter P{^) or also the famous Newton's problem of the body of minimal resistance 
as studied by T. Lachand-Robert and coauthors: see for example ||5l [161 and see also |[T4l |71 for more examples 
arising in the operator theory. Actually, the techniques employed in |[T4l . and here as well for (ii), are inspired 
from those introduced in llT6l . The main novelty here in the results (ii) is that the functional are not necessarily of 
the local form Q and may include shape functional defined through state functions which are solutions of partial 
differential equations (PDE). The "concavity condition" is then expressed in a functional way through the coercivity 
of the second derivative in an adequate functional space : see Theorem l2.9[ In 121, a similar concavity phenomenon 
is used to get qualitative properties of minimizers in higher dimension, under assumptions about their regularity and 
convexity. We avoid here any assumption of this kind for the planar case. 

The general optimality conditions including the infinite dimensional Lagrange multipliers were also provided 
(and exploited) in the same paper lfT4l . They are revisited here in an VF^'°^ -context which is better adapted to our 
more general functional (see e.g. Proposition [ST]). 

Similar arguments to those used here to obtain the first family of results (i) may also be found in |l3| where op- 
timality conditions with convexity constraints are developed in an A^-dimensional setting. They are exploited for 
several examples in dimension 1 (or in radial situations) to obtain C^-regularity of the optimal shapes. With our 
approach here, we are able to reach -regularity and this is valid for a rather general family of functionals. 

About a localization of the approach: Let us mention that our two families of results may be mixed in the same 
functional: indeed, as often the case, it may be that the required convexity property for (i) is valid on some part of the 
boundary of the optimal shape, while the concavity property for (ii) is valid on the other part. Then, the techniques 
developed here may be locally applied to each part and we can obtain at the same time smooth and polygonal pieces 
in the boundary. However, as one expects, it remains difficult to understand the portion of the boundary which 
remains at the intersection of these two parts. We refer to Section lSTTl for more details. 

To end this introduction, let us say that many questions are of interest in shape optimization among convex sets. 
Here, we try to exploit as much as possible analytical tools to obtain precise qualitative results for optimal shapes 
among convex planar sets. But many questions are left open in higher dimensions. Among them, and besides the 
Newton's problem already mentioned, we can quote the famous Mahler conjecture about the minimization of the 
so-called Mahler-product |i^||i^°| among symmetric convex bodies in M'^ (see llT9l ). which is of great interest in 
convex geometry and functional analysis, and the Polya-Szego conjecture about the minimization of the Newtonian 
capacity among convex bodies of R'^ whose surface area is given (see for example H and reference therein). 

This paper is structured as follows. In the following section we state our main results. In Section [3] we focus on 
the regularity result (i) and we apply it to some various examples. In Section IH we deal with problems leading to 




(6) 
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polygonal solutions (result (ii)), and we again consider in detail some classical examples. We conclude with some 
remarks and perspectives. 

2 Main results 

2.1 Notations and problems 

We set T := [0, 27r). Throughout the paper, any function defined on T is considered as the restriction to T of a 27r- 
periodic function on M. We define W^'°^{T) := {u G wl^^{M.),u is 27r-periodic}, and similarly for any functional 
space. If u G VF^'°°(T), we say that u" + u>Oif 

\/ve W^'^^iT) with v>0, I {uv- u'v') dO > 0. 

In this case, u" + u is a nonnegative 27r-periodic measure on M and finite on [0, 27r]. 

We denote by Sad a class of open bounded sets in M? (including constraints besides convexity). We will focus on 
two problems: 

min{J(r2), Q G Sad, biconvex}, (7) 

min{J($7), ft £ Sad, ^ convex, M(r2) = Mq}, (8) 

where J : Sad — K is referred as the energy and M : Sad — ^ is an extra constraint (Mq given in G M"'). 

In order to analyze the regularity of an optimal shape, we transform these problems into minimization problems 
in a functional analytic setting as follows: choosing an origin O and using parameterization ([2]), we define 

Tad-={ueW^'°°{T), n^eSad}, (9) 

the set of admissible gauge functions, endowed with the || • ||^yi,oo(j)-norm, and we assume that this set can be 
written 

j;^ = {n G Ty^'~(T) / A;i < u < /C2 and u > 0}, (10) 

for some functions ki,k2 : T — M+ respectively upper- and lower-semicontinuous (see Remark [2?T] below for this 
assumption). 

A simple calculus of the curvature shows that is convex if and only if u" + u > 0. Moreover, the support of 
the measure u" + u gives a parametrization of the "strictly convex part" of the boundary, and a Dirac mass in this 
measure correspond to a corner of the associated shape; we have for instance that is a convex polygon if and 
only if u" + u is a finite sum of positive Dirac masses. 

If is a solution of problem ^ (resp. ([8])), then its gauge function uq is respectively solution of: 

j{uo) = inm{j{u), u" + u>0, Tad}, (H) 
resp. j(^o) = min {j{u) / u G J-'ad, u" + u > and m{u) = Mq}, (12) 
where j : Fad ^ := J{Qu), and m : Tad — ^ 1^, m(n) = M{Qu)- 

Our main goal in this paper is the analysis of the convexity constraint. Thus, given an optimal shape Qq, we focus 
on the part of dQo which does not saturate the other constraints defined by Sad- We therefore define, for uq G Tad 
and ^0 = 

Tin ■■= TiniTad,Uo) = {e £T / ki{e) <Uo{e) <k2{e)}, (13) 
{^no)^n := \xedno/^OeTin,X = ^— {cose, sine)]. (14) 

See Example |2!2l and Figure 1 for examples. 
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Remark 2.1 If ki or k2 happened not to be semicontinuous, we could replace them by 

ki = mi{k : T — )• M continuous ,k>ki}, /c2 = sup{/c : T — )• R continuous ,k<k2} 

and we have 

{u G W^'°°{T) / ki<u<k2} = {ue VF^'°°(T) /h<u<k2}- 

Therefore, the assumptions on ki and ^2 are not restrictive. Note that, thanks to the regularity of uo,ki,k2, the set 
Tin is open. 




Figure 1 : Inclusion constraints 

Example 2.2 A frequent example for admissible shapes Sad is: 

Sad ■= bounded open set of / K2 C C i^i } , 

where K2 and Ki are two given bounded open sets. If for example Ki and K2 are starshaped with respect to a 
common point O, chosen as the origin, then 

J'ad = {ue W^'°°{T) /ki<U< k2}, 

where ki, k2 are the gauge functions of K2 and Ki respectively. In that case, given a set G Sad, 

{dn)in = dn\idKi\jdK2), 

see Figure 1 . 

The analysis of the optimal shape around the set {0 / uo{0) G {ki{9), k2{0)}} = T \ Tj„, where the inclusion 
constraint is saturated, may require more efforts, see |[T4l for example. In this paper, we will not discuss this 
question. 

Note that we can also consider the case K2 = and/or Ki = with ki = and ^2 = +00. 
Example 2.3 With respect to the constraints m, M in (fT2]) . a classical example is the area constraint: 

m(n) := = ^0 ^^=^ j = ^0, 

where \Q\ denotes the area of O. 
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2.2 The main results 



As explained in the introduction, Section [T] we will prove two types of results: they are described in the two 
following subsections. 

2.2.1 "Optimal shapes are smooth" 

First we consider the problem dV]) and its associated analytical version (fTTT) . We assume that J{Q,) = R{Q,) + C(il), 
R satisfying some "regularity" assumption, and C being written like in and satisfying a convexity like property. 
More precisely: 

Theorem 2.4 Let uq > be an optimal solution of (1111) with J^^d of the form dlOl ) and 



where r and G satisfy: 

i) r : tyi'°°(T) M is around uq and G : {9, u,q) e T x {0,oo) x M. R is around T x no(T) x 
Conv(ug(T)), where Conv(Mo(T)) is the smallest (bounded) closed interval containing the values of the 
right- and left-derivatives Uq{6~^),Uq{6~),6 G T, 

ii) r'{uo) £ LP (T) for some p G [l,oo], 
Hi) Ggg >OinT X uo{T) x Conv(u^(T)). 



See Section [3?T] for the proof, and Section[3]2]for explicit examples. 

Remark 2.5 A -regularity result has been proved for a similar problem with r = in |[3l with different boundary 
conditions, with a proof which is also based on first order optimality conditions. Here, for periodic boundary 
conditions (but this is not essential), we improve this result to the C^'^-regularity, and generalize it to the case of 
non-trivial r, which is of great interest for our applications. Let us also refer to [6] for a higher dimensional result. 

Let us remark that the same result is valid, with the same proof, if we only assume that r'{uo) is the sum of a 
function in ^^(T) and of a nonpositive measure on T. □ 

We can also get a similar result for the equality constrained problem ^ and the associated problem (fT2l ) as follows. 

Theorem 2.6 Let uq > be an optimal solution of (1121 ) with j, J^ad in Theorem \2.4\ and m : W^'°° — ?• M'^ a 
function around uq with ■m'{uo) G (L'P(T))"' onto. Then 




(15) 



Then 



Uq G W' 



^'^(Ti„), where Tin defined in (1131 ). 



Mo e W' 



{Tin). 



See Section [3?T] for the proof, and Section[3]2]for explicit examples. 



For a shape functional, using parametrization Theorems 12.41 and |Z61 lead to the following. 
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Corollary 2.7 Let Sad be a class of open sets in such that Fad '■= {u / Q.u ^ Sad^ is of the form (llOb (Q.^ is 
defined in Q ), and let J : Sad -^"^be a shape functional: 

i) Let Oq be an optimal shape for problem and assume that J = R + C with: 

Vu G J^ad, R{^u) = r{u) and C(17„) = / 0(9, u{e),u'{e))de, 

Jr 

where r and G satisfy assumptions ofTheorem \2.4\ for some p G [1, oo]. Then {dQo)in, as defined in (|14l) . is and 
its curvature is in L*'((9r2o)m)- 

ii) A similar results holds for the problem ifm{u) = M{Vtu) satisfies the hypotheses in Theorem \2.6\ 

Remark 2.8 The results of this section are in an abstract analytical context, and do not depend on the charac- 
terization of the domain. Therefore, one could consider the classical characterization of a convex body with its 
support function instead of the gauge function. In Section [l!4l we give a geometrical interpretation of similar- results 
associated to this parametrization. 

2.2.2 "Optimal shapes are polygons" 

Our second result is a generalization of Theorem 2.1 from |[T4l . We give a sufficient condition on the shape func- 
tional J so that any solution of ([T]) be a polygon. In llT4l . the first two authors only consider shape functionals 
of local type like The following results deal with non-local functionals, which allow a much larger class of 
applications, including shape functionals depending on a PDE. 

Theorem 2.9 Let uq > be a solution for (fTTI ) with Fad of the form (ITOl) . and assume that j : 1^^'°°(T) — )• M is 
around uq and satisfies (see Section \4~l\ for definitions of H^-(semi-)norms): 

3s G [0, 1), a > 0, /3, 7 G [0, oo), such that 

yv G W^'°°(T), j'\uo){v,v) < -a|?;|Hi(T) +7l^l//i(T)lbk-(T) +/3||v|Ih=(t)- (16) 
If I is a connected component of Tin (defined in (1131 ) ). then 

Uq + no is a finite sum ofDirac masses in I. 
See Section l4~T] for a proof and Section 1421 for explicit examples. 

Remark 2.10 We can even get an estimate of the number of Dirac masses in terms of a, /3, 7, see Remark 14^21 

Remark 2.11 Theorem 12.91 remains true if (fT6] ) holds only for any v such that (denoting fi = Uq + uq): 

3ipe L°°(T, p.) with v" + v = Lpfi. 

Indeed, the proof of Theorem 12.91 uses only this kind of perturbations v which preserve the convexity of the shape. 

□ 

As in Section [2.2. II we can also handle the problem with an equality constraint as follows. 

Theorem 2.12 Let uq > be any optimal solution of (|12l) with j, Fad cis in Theorem \2.9\ and the new assumptions: 
j'{uo) G (C°(T))' , and m : W^'"^ R'^ is around uq, 

m'(no) G {C^{T)Y is onto, \\m" {uq){v,v)\\ < P'\\v\\]js(jy for some p' G M. 
Then, if I is a connected component of Tin (defined in (1131 ) ). 

^0 + ^0 is a finite sum ofDirac masses in I. 
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See Section l4T] foi' the proof. 



Again, using the parametrization we get the following result. 

Corollary 2.13 Let Sad be a class of open sets in such that Fad '■= {u / Vtu ^ Sad} of the form (1101 ). l^o 
an optimal shape for the problem ©for ©/or the constrained problem), and assume that j : u £ J- ad ^ J{^ii) 
satisfies assumptions of Theorem 12. 91 ( and m : u £ J-ad ^ ^{0,^) satisfies assumption in Theorem \2. 12\ in the case 
of the constrained problem). Then: 



Remark 2.14 When one uses the parametrization of convex sets by the gauge function u, Vlu is a polygon if and 
only if u" + u is a sum of Dirac masses. When parametrizing Q. with the support function as in Section [331 one has 
the same characterization. Therefore, the results of this section hold if we work with the optimization problems as 
in Section [Ml 

3 Shape functionals containing a local-convex term 

In this section, we give the proof of the results in Section [2.2.11 that is to say regularity results for solutions of ([TTI) 
or ([T2I ). Using the parametrization ©, since the regularity of a shape and of its gauge functions are the same, we 
consider several applications of regularity for optimal shapes to classical examples of energies. We conclude with a 
few remarks about the application of our results when we use another parametrization of convex bodies, namely the 
support function. In that case, we get the regularity of the support function, which does not imply the regularity of 
the corresponding shape, but only the fact that this one is strictly convex. 

3.1 Proof of Theorem [U and [Ml 
First order optimality condition: 

A first optimality condition for the problem ([TTT l is stated in |[T4l Proposition 3.1, 3.2] when j is defined and differ- 
entiable in the Sobolev Hilbert space ff^(T). We give here an adaptation to state this result in instead (which 

is important for our applications involving a PDE, since the shape functionals aie known to be differentiable for 
Lipschitz deformations only). 

Proposition 3.1 Let uq > be a solution of f[771) with j : VK^'°°(T) — )• M of class and such that j'{uo) G 
C°(T)'. Then there exists Co e VF^'^fT), such that 



Remark 3.2 Without any assumption on j'{uQ), we would a priori get a Lagrange multiplier Co S L°°(T) (see the 
proof below). The non-continuity of Co may lead to some difficulties, especially to state that Co = on Suppfug + 
Uq). Though a restriction, the assumption j'{uo) G C*'(T)' will be satisfied in all of our applications. 

Proof. We set 



each connected component of (9r2o)m polygonal. 



Co > onT, Co = on Supp(tiQ + uq), and 



< 




(17) 



g:ve W^'°° ^v" + ve {W^'°°y in the sense that {v" + v, if) 




8 



and we consider Y := Im(g) = {/ G W^'°°{Ty, (/, cos)(iyi,°o)/xVKi.°° = if: sm)(^^n,ooy^^rl,cx, = 0}, which is a 
closed subspace of (VF^'°°(T))'. 

Applying the same strategy as in |[T4ll . one gets Zq £ such that lo{g{uQ)) = and 

V/ey, />0^/o(/) >0, and Vt; G (j'K), 7;)(H.i,ooy = (^o, + ^^>y'xy ■ 

We restrict ourselves to G V{T) := C°°{T), and consider 

Co : / G i^(T) n y ^ (Co, fh'xv ■■= ik, Dy'xy- 

Our aim is to prove that Co can be extended to a continuous linear form on L^(T). First, for / G D(T) n y = {/ G 
P(T), Jrjp/sin = f rp f cos = 0} we choose the unique f G W^'^{T) such that {Jrjpusin = J^u cos = 0} and 
v" + V = f in T. Then there exists C < oo independant of or / such that 

lbllvKi.°°(T) < C||/IIli(T)- (18) 

Indeed, we first get an L°°-estimate using Fourier series: if / = J2n& f{''^)^n with e„(0) = e*"^ and f{n) = 
/(^)e-^"^f , then v = T3^i^/(n)e„, and therefore 

max|/(n)|<C||/||ii, 




with C < oo. Then we get a VF^'°° -estimate by choosing 9q such that v'{Oq) = (which is always possible, thanks 
to regularity and periodicity of v), and getting from v" + v = f that 



{f{s)-v{s))ds 



<27T{\\v\\L^ + \\f\\Li 



\v\e)\ - 

which concludes the proof of the estimate (fT8l ). 

Therefore, we can write (C may define different universal constants) 

V/ G y nP(T), KCo,/>I''xd| = \{l,v" + v)y'xy\ = \{j'{uo),v)(^w^,^YxW^^^\ < C\\v\\wi,^ < C||/||ii. 

(19) 

We now extend Co on V{T) by 

V/ G V{T), (Co, f)v'xv = (Co, / - 7(l)ei - f{-l)e-i)v'xv. 

Then, applying ^ to / - /(l)ei - /(-l)e_i, we get 

V/ G P(T), KCo, fh'xvl < C\\f - /(l)ei - /(-l)e„i||ii < CWfhi, 

and therefore by density, we extend Co to a continuous linear form in L^, which can be identified with Co £ 
Moreover, in the sense of distributions: 

{Co,v" + v)v'xv = {j'{uo),v)v'xV, that is to say Co + Co = j'{uo)- 

From the hypothesis for /(no) it follows Co + Co £ (C'^(T))' which implies Co G W^'°°{T). Using the continuity 
of Co and the fact j'{uo){uo) = we get fj Cod{uQ + no) = by a density argument. Therefore, the rest of the 
proof stays as in |[T4l . namely we prove that we can add a combination of cos and sin to Co so that Co ^ 0. □ 
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Proof of Theorem |231 

Applying the previous proposition, and using the hypotheses on the functional j, we get: 

yv£C^{Tin), f{uo)v = r'{u)v+ Gu{6,uo,UQ)v + Gg{e,uo,UQ)v' = (Co + Co ; ^)(lyl><-{T))'xM/l><-(T)■ 
JT 

To integrate by part in this formula, since u'q is only in BV{T), we may look in [201 (see also HI) to get: 

r'{uo) + Gu{0,uo,Uq) - Geg{0,uo,UQ) - Guq{0,uo,u'Q)uQ - ulGqq{6,UQ,UQ) = Co + Co inP'(Ti„). (20) 

where Gqq{6, uq, u'q) = Jq Gqq{0, uq{6), (1 — t)tiQ(0+) + tu'Q{9~))dt. For simplicity, we will drop the indication 
of the dependence in (0, uq,u'q) and write more simply 

r'(no) + Gu- Geq - Guqu'^ - u'^fi^q = Co + Co in V'{Tin). (21) 

Equality (|2TI ) implies that Co is a Radon measure, and also that the singular parts of the measures in the two sides 
of (I2TI ) are equal. To study the sign of these measures, we will use the following lemma. 



Lemma 3.3 The measure Co satisfies: Co > on [Co = 0]. 

Proof of Lemma|331 Let G (M) , v? > and let p„ : M+ M+ be defined by 

Vr G [0, 1/n], Pnir) = 1 — nr; Vr G +cx3), Pnir) = 0. 
Recall that Co G ^^^'^(T) and Co > 0. Then 



/ 



Wn^diO = - (1 C^9^Vn(Co) +V=K(Co)C^') >- j C'o^'PniCo) 



Letting n tend to +00 leads to 



vdiO > - / Co<^' = 0, 

[Co=0] ^[Co=0] 



the last integral being equal to thanks to the known property Co = a-^- on [Co = 0]. □ 
End of the proof of Theorem g^l 

Denote K := Supp(no + uq). Recall that Co = on K by Proposition [37T] By Lemma [331 Co > on K. Let 
= fJ-ac + f^s and Co = ?^ac + bc the Radon-Nikodym decompositions of the measures Uq , Cq in their absolutely 
continuous and singular parts. Note that: [uq + uq > ^ Hs > 0] and > on K. 

Identifying the singular parts in the identity (|2TI ). and using that r'(uo), Gu, Ggq, GuqUQ, u^Gqq are at least L^- 
functions, we are led to —j-LgGqq = Ug in Tj„. Since Gqq > 0,fis > 0, > on D Supp(/_is), we deduce 
Us = = Ug in Tin- Thus, uq G VF^'^(Tj„) and u'q is absolutely continuous on Tj„. In particular, Gqq = Gqq on 

in- 

We can now obtain higher regularity, using again the multiplier Co ■ Indeed, on one hand, we deduce from Lemma 
13. 3[ from (|2T]) and from the inequality — UgGgq < u^Gqq, that, on the set Tj„ n K 

< Co < '^'(^O) + Gu - Geq - GuqU'Q + UoGqq G 1^(7). 

Thus, Co £ L^O^in n K). Going back to (|2TI ) and using that Gqq = Gqq is bounded from below on the compact set 

T X uo{T) X Conv{u'Q{T)), we deduce u'^ G LP{Tin n K). 

On the other hand, in the open set Tj„ \ K, we have u'q + no = so that u'q G L°°(Tj„ \ K). As a conclusion 

< G LP{Tin). □ 
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Proof of Theorem IIS 

Optimality conditions are written with the Lagrangian (since m'(no) is onto, see also lfT4l Proposition 2.3.3]): 

\/v G C^{Tin), j'{uo)v + fi- {m'{uo)v) = (Co + Co ' ^^>(vyi.-')'xW'i.°°, 

for some G M"'. The regularity of m'{uo) implies that the strategy used in the proof of Theorem I2.4l remains valid. 

□ 

3.2 Examples 

In this section, we apply Corollary 12.71 to a number of classical energy functional. For the proof of the differentia- 
bility of the shape functionals see Section [331 We start by reminding some classical PDE functionals that we use in 
our examples. 

Dirichlet energy - Torsional rigidity 

For an open bounded set in M^, we consider the solution of the following PDE, in a variational sense: 

UneH^iQ), -AUn = finQ, (22) 
and we define the Dirichlet energy of O by 

About the regularity of the state function, we are going to use the following classical result (see |[T3l . |[9l). 

Lemma 3.4 Let Q, be convex, f £ Lf^^O^"^) with p> 2, and Un be the solution of^. Then Un G W^'^^iQ) n 

Remark 3.5 When / = 1, the Dirichlet energy is linked to the so-called torsional rigidity T{Q), with the formula 

T{n) = -2Ei{n). 

First Dirichlet-eigenvalue of the Laplace operator 

We define Ai(r2) as the first eigenvalue for the Laplacian with Dirichlet's boundary conditions on d^. It is well- 
known that, if we define Uqhs a. solution of the following minimization problem, 

Ai(0):= / \VUn\^=nim{ [ \VU\^ U e H^{n), [ = l] , 
Jn Uq Jq J 

then Un is (up to the sign) the positive first eigenfunction of — A in fi: 

Un G H^{n), -AUn = Xi^Un, [ U^ = l- 

Jn 

Again, like in Lemma[331 if n is convex then Un G H'^{n) D W^'°°{n) and C/f^ > in n. 
We are now in position to state some applications of Corollary 12.71 
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Example 3.6 (Penalization by perimeter) One can study 

mm{J(0) := F Ef{n), + P{n) / n convex, Di C n C D2} (23) 

where F : (0, +00) x (-00, 0) x (0, 00) R is C\ / G Di,D2 are bounded open sets, Ef{Q) is the 

Dirichlet energy and Ai(0) is the first eigenvalue of — A defined as above. 

Proposition 3.7 I/Qq is an optimal set for the problem (|23] ). then the free boundary dQo n {D2 \ Di) is C^'^ (or 
equivalently W'^'°^), that is to say dO.^ n {D2 \ Di) has a bounded curvature. 

The proof is a simple consequence of Section [331 which asserts that R{Q) = F{\Q.\, Ef{Q), Ai(r2)) and C{Q,) = 
P(0) satisfy the assumptions in Corollarv 12.71 with y = 00. 

Note that in Proposition 13. 71 we could also add a dependence of F in the capacity of $7 or in any shape functional 
which is shape differentiable and whose shape derivative can be represented as a function of L°°{dQ.) when Q, is 
convex. 

Remark 3.8 The constraints Di C ^ C D2 helps existence for the problem (1231 ). Of course, if one can prove 
existence of an optimal shape without these constraints (mainly, one need to prove that a minimizing sequence 
remains bounded and does not converge to a segment), the result of Proposition I3.7l remains a fortiori true for the 
whole boundary of the optimal shape, i.e. OUq is C^'^. 

Example 3.9 (Volume constraint and Perimeter penalization) We can also consider a similar problem with a 
volume constraint: 

mm{J{n) := F{Ef{n),Xi{n)) + P{n) / convex, and \n\ = Vq}, Vq G (0,+oo). 

In this case, the first optimaUty condition will be similar to the one for the problem (1231 ) with F{Ef{Q), Xi{0,)) + /u|r2| + P{0,), 
where /Li is a Lagrange multiplier for the constraint = Vq. Theorem 12 . 6 1 applies and one gets globally the same 
regularity result (but global) as in Proposition 13. 7l on any optimal shape. 

Example 3.10 (Perimeter constraint) If one considers again a problem with a perimeter constraint, 

mm{J{n) := F {\n\, Ef{n), Xi{n)) / n convex, and = Pq} (24) 

where Pq G (0, +00), one needs to be more careful. In this case, the first optimality condition will be similar to the 
one for the problem (|23] ). with F{\Q\, Ef{Q), Ai(r2)) + p,P{Q,), where is a Lagrange multiplier for the constraint 
P(ri) = Pq. Therefore if we are able to prove jj, > then we can apply the same strategy as in Theorem 12.41 and 
we therefore get the same regularity result as in Proposition 13.71 However, if ^ < 0, we refer to Example 14.91 

Example 3.11 In a more abstract context, one can consider 

min{ J(0) - a\n\ + P{n) / n convex C D}, (25) 

where J is a shape differentiable functional, increasing with respect to the domain inclusion, D is an open set, and 
a > (if a = 0, the empty set is clearly solution of the problem). Again, we get that d^^o n D has a locally bounded 
curvature. Indeed, the derivative of j{u) := J(ilu) is a nonpositive measure, thanks to the monotonicity of J (see 
lITSl ). and we apply Theorem |23] combined with the end of Remark |23] 

3.3 Computation and estimate of first order shape derivatives 

In this section we will prove the differentiability of the shape functional involved in the examples of Section [l!2l 
which aie needed in Proposition 13.71 
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3.3.1 Volume and perimeter 

About geometrical functionals, it is easy to write the area and the perimeter as functional of u, namely 



a(n) := = / ^d9, p{u) := = / ^^^^^d9, u G W^'^{T) r\{u> 0}. (26) 

I 2 I 2 

Note that p(n) = JjG{6,u{9),u'{9))d6 withG{6,u,q) = ^ and one can easily check that Gqq = ^^2_)_g2)3/2 > 
0. " " ^ 

3.3.2 Dirichlet Energy - Torsional rigidity 

We focus our analysis around a convex open set JIq with parametrization uq > 0. For \\u — uo\\-^^i,ac(j^ small, 
consider 

e/ : W^'°°{T)n{u>0} R, 

In order to study the differentiability of ej near uq, we use the classical framework of shape derivatives. As 
usual, we need to work with an extension operator: the deformation d^Q to d^u allows to define the vector field 
^(u) : — ^ ^ such that d^u = {Id + ^(u))(5r2o)- We will consider an extension to of this transformation, 
since we need to study the differentiability of u — )• := Uq_^^ o [Id + i{u)) G i7g(J7o)> where Uu '■= Uq^ (see 
ifTOl for example). 

If we consider a smooth extension operator ^ : W^'°°{T) W^''^{M.'^; M^), we have {Id + ^{u)){dno) = dn^ 

if 

where {:^{9),9) are polar coordinates (for simplicity, we will often write uq, m or ^ instead of uq{9), u{9) or 

m{r,o)). 

Remark 3.12 The transformation ^(u) can be extended to R? in different ways. The easiest way is to take 

i{n){r, 9) = (J- - e'%{r, 9) in M^, (28) 

\u{9) uq{9)J 

where 77 G C^(M^), = in a neighborhood of the origin and = 1 in a neighborhood of dQ.o. 

This (polar) extension of i{u) is such that ^ G (;°°(W^^'°°(T); W^1'°°(M2; ^2^) near mq, and is sufficient for the 
results of this section. More work will be needed for the second order shape derivatives, see Section l4.3.2[ 

Let us point out that if ^ is in a neighborhood of uq and satisfies (l27l) . then 

2 

yveW^^^{T): i'{uo){v) = -^e''^, i" {uo){v,v) = 2%e'^ on 8^0- (29) 

Note also that the method used in the proof of Lemma D.14[ which is needed in the proof of Proposition |3.13[ allows 
to say that the method a priori fails if we consider an extension operator ^ : II^{T) — >• H^{M?; M^). This explains 
our choice to work with v G W^''^{T) rather than v G H^{T), even though it introduces extra difficulties (like in 
Proposition 13.11 and in the proof of Proposition 14.111 ). □ 

The main result of this section is the following. 



13 



Proposition 3.13 Let^Q = convex, f £ Hj'^^iM.^), k e W and £, e C''{W'^'°°{T);W'^'-°°{M.^;M.'^)) near uq. 
We have: 

i) e f is near uq. 

ii) If satisfies (l27l l, then for any v G W^'°°{T) we have 

e'j{uo){v) = - [ hvUo\He{no){v)-i^o)dso= [ hvUo{xe)f47lde, (30) 
whereUo € (^n) is the solution of (1221 ) in Qn. fn is the exterior unit normal vector on d^n, xg = mj(g) (cos 0, sm.9) € 
Hi) Furthermore, ej(tio) G L°°{T). 

The proof of this proposition is classical and uses the following lemma, which will be needed in the following 
section. 

Lemma 3.14 Let uq G W^^°°{T), uq > 0, f € i^L(I^^)' ^ ^ I^*- We have: 

i) The map u G W^'°°{T) ^ Uu e F^f^o) is near uq. 

ii) Forve W^^°°{T), set 

UL{uo){v), C/^ := U;,-VUo-^'{uo){v). (31) 

Then 

U'q(^L^{^q), AC7^ = mP'(Oo), (32) 
% + VU^-i'{uo){v) G Hl{^^). (33) 

Furthermore, if Uq + > 0, then Uq G H^{Qq). 

Remark 3.15 Here we are not interested in the differentiability of u i— >■ C/„ and the function Uq is directly defined 
by (OTI ). In fact, the map u ^ Uu (with C/„ extended by zero in M^) is differentiable in L^(M^) and its derivative 
equals Uq in Q.q, see Theoreme 5.3.1, lITOl for example. 

Proof of Lemma l3.14t 

i) The map 9 G ty^'°°(R2; M^) ^ U^ja+e){no) ° i^d + 9) e H^{no) is C'' in a neighborhood of 0, see for 
example ifTOl Proposition 5.3.7]. We conclude by using the composition of this map with ^. 

ii) It is clear- that U(^ G L^i^o) and that U(^ + VC/q • £,'iuo)iv) = Uq £ H^i^o)- To prove A?/^ = we 
consider the map S : VFi'~(T) ^ W'^''^(R'^;R'^), S{u) = {Id + ^{u)^, which is well defined and C'' in 
a neighborhood of no. From S{u) o {Id + £,{u)) = Id, it is easy to check that for v G W^'°°{T) we have 

S'{uo)iiv) = -('iuQ)iv), S"iuQ)iv,v) = 2V^'iuQ)iv) ■ C'iuQ)iv) - ^"iuQ){v,v). (34) 

Let ip G V{Qq). From ((22]) . for all u near uq we have J^^^ 11^ o S(u)Aip — fip = 0. Differentiating this 
equality on the direction v gives 

j (U^ o Siu) + VUu o S{u) ■ S'iuo)iv)^ A^ = 0. (35) 

Replacing n = in (1351 ) and using (l34l ) gives 



which proves ii). 
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iii) If Uq + uq > then is convex. From Lemma we obtain Uq G H'^{Q.q), which implies U'q G H^{Q.q). 
□ 



Proof of Proposition 13. 13 1 

i) The functional u i— >• ef{u) can be seen as ej(ti) = <S'/,no ° ?(^)' where <S'/,no is a classical functional, introduced 
to compute shape derivatives: 

e ^ Ef{{id + e){^o)). ^ ' 

As ^ is near uq and £:/,f7(, is near = in Ty^'°°(R2; M^)^ see US Corollaire 5.3.8]), the differentiability of 
ef{u) follows. 

ii) As we have ey(u) = — ^ Uu ° S{u)f and C/„ = on (?r2o> from Corollaire 5.2.5, jlOI . we obtain 

e'j{u){v) = (u',oSiu) + VUuoSiu)-S'iu)iv)') f. (37) 
Taking u = uoin the last equality and using (l34l ) gives 



Finally, by changing the variable sq = 'v/"o+J"o) _^g^ taking into account that z/q = (— e*^ + ^(ie*^)^ — ^ 



and after using ( |29l ). we obtain (l30l ). 

iii) As G iV* it follows / G LP{no), for all p G [l,oo). Then Lemma [H gives Uq G W^'°°{no), so e^(no) G 
L°^(T). □ 

3.3.3 First eigenvalue of the Laplace operator with Dirichlet boundary conditions 

We consider 

li: {u G W^'°°{T), u > 0} ^ M 

n I— )• Zi(n) := Ai(r2u) 

and we have the same result as in Proposition 13. 131 see for example Theoreme 5.7. 1, ifTOl . and (l29l) . with 

Jt W) 
3.4 Application with the dual parametrization 

Instead of using parametrization by the gauge function, one can also use the well-known parametrization by the 
support function of a body, namely 

ye G T, hn{9) := max{x • e'^ , x G n}. 

We get a characterization of the convexity in a similar way to (O: 

Q, is convex =^ + Hq > 0. 

Conversely, if /i G W^'°°{T) satisfies h" + h > 0, then one can find a unique (after a choice of an origin) open con- 
vex set, denoted 0^', whose support function is h (see lITSll for example). This parametrization is the dual of the one 
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with the gauge function. Indeed, the gauge function of Q is the support function of the dual body of and vice versa. 



Therefore the optimization problem 

min{ J(il) / $7 G Sad, ^ convex}, (38) 
where Sad is a class of open planar sets, becomes 

find ho G Tad such that j(/io) = niin{j(/i), h G Tad, h" + h> 0}, where 

m = J{n^), and Kd = {he W^'^iT) / G Sad}, 
which is the same as (fTTI ). 



(39) 



Again, if the set of admissible functions can be written 

T^d = {he l^i'°°(T) /h<h< k2}, (40) 

we can define = {6 e T / ki{e) < h{e) < /c2((9)}, and then (30)^ = {x e dn s.t. 39 ^%a, x • e*^ = h(Q)\, 
i.e. the set of points of 50 whose supporting plane is orthogonal to (cos(0), sin(0)) with B G Tj„. 

As in Example 12.21 if Sad = {O / A'l C O C K2\, where K\ and are two convex open sets, then (l40l ) is 

satisfied with fci, k2 the supports functions of A'l, K2, and in that case (90)^ = <90 \ (5/^1 U dK-}). 
Therefore one gets a dual version of Corollary I2.7l as follows. 

Corollary 3.16 Let Oq = be an optimal shape for the problem (1381 ) with J = R + C, and assume that, 

yheT^d, R{^^) = r{h)andC{^t)= [ G{9,h{e),h'{9))de 

where r and G satisfy the assumptions o f Theorem \2.4\ for some p G [1, 00]. Then 

ho G W^'PifTn). 
This implies in particular that ((90o)in strictly convex. 

Remark 3.17 This parametrization is especially interesting when one has to deal with the perimeter because in this 
case P{^}^) = Jjhd9. An example of a function C{fi^) satisfying the hypotheses of Corollary 13. 161 is now the 
opposite of the area, since 

2 Jt 

However, it is not easy to work now with functionals coming from PDE. Indeed, it is well-known for example, 
that the derivative of Ai in terms of h is not more regular than a measure on T, see |[TTl[T2l . We think that this can 
be explained by the fact that some solutions of problems like (l23l) may not be strictly convex. 



4 Optimization of concave non-local shape functionals 

In this section, we prove the results of Section 12.2.21 The main proof relies on the analysis of the second order 
shape derivatives. Next we apply these results to various energy functionals involving the Dirichlet energy or the 
first eigenvalue of the Laplace-Dirichlet operator. Since the optimal shapes come with no a priori regularity except 
the convexity condition, one needs some delicate computations to check the required assumptions. This leads to 
rather sharp estimates on second derivatives which are interesting for themselves. 
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4.1 Proof of Theorems US] and HH 

We first introduce tiie classical Sobolev semi-norms on T. For s G we set: 

WlH^m := where u(n) := [ ti(6')e"*"^ — . 

ngZ 

We also define H'^{T) := {u G L^(T) such that \u\fj3(j-^ < +00} and •= ll^lli^m + 1^ 



Proof of Theorem IIS 

The main idea is to prove that for a deformation supported by a small set, the estimate ([T6l ) is a concavity estimate, 
and so it violates the second order optimality condition. This relies of the following Poincare-type inequality: 

Lemma 4.1 Let s G [0, 1) and e G (0, vr). Then there exists a constant C = C{s) independant on e such that, 

\/u G H^{T) such that Su]ip{u) C [0,e], ||?i||//ii(T) ^ Ce^~^\u\jji(jy 

Proof of Lemma \4A\ Let u G C°°(T) with Supp(n) C [0, e]. If we first assume that s = 0, then we have the 
classical Poincare inequality (with the optimal constant), proved using the fact that Itil^i^-j.^ = Jj u'"^, so 

If one has now s G (0, 1), one can proceed with an interpolation inequality, easily obtained by Holder inequality: 



and so 



l-s 

= |np1n(n)p«|n(n)p{i-) < ( ^ |n|2|n(n)|2 ) ( ^ |n(n)'2 



□ 



Let K := Supp(nQ + uq). Assume that, for a connected component / of Tj„, K n / is infinite. Then, there 
exists 6*0 G I an accumulation point of ii' n /. Without loss of generality we can assume ^0 = and also that there 
exists a decreasing sequence (e„) tending to such that i^' n (0,e„) C / is infinite. Then, we follow an idea of 
T. Lachand-Robert and M.A. Peletier as in lfT4l (see also |[T6l ). We can always find < < e„, « = 1, . . . , 4, 
increasing with respect to i, such that Supp(Mo + ^^o) H (e^, e^^'^) / 0, i = 1, 3. We consider VnA G Vl^^'°°(T) 
solving 

Vn,i + '"n,i = l(^j^^gj+i)(no + no), Vn,i = in {Q,enf, i = 1, . . . , 3. 
Such Vn,i exist since we avoid the spectrum of the Laplace operator with Dirichlet boundary conditions. Next, we 
look for An,i, i = 1, 3 such that Vn = ^ Xn,iVn,i satisfies 

i=l,3 

The above derivatives exist since Vn^i are regular near and e„ in (0, We can always find such A„.i so as they 
satisfy two linear equations. It implies that v'^ does not have any Dirac mass at and e„. It even implies that the 
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support of Vn is included in [e^, e^]. In particular, v" + f.„ = ip{uQ + uq) where is bounded and with support in 
[e\, e^j. As Supp(no) n {el^,el^^) ^ 0, we also have v„ / 0. 

Since Supp(f„) C Tj„ and f + i^n = fi^o + ^*o), it follows that no + tvn is admissible for |t| small enough 
(and n fixed). Consequently, since j{uQ + tVn) > j{uo) for \t\ small, we have j'{uo){v) = and then by using the 
assumption ([T6l ) and Lemma l4~n we get 

< j"{uo){Vn,Vn) < -a| I'n ll^i (T) + T^n | (T) l^n (T) + I^n || (T) (^1) 

< (-a + C77e^-^ + C2/3(e„)2(i-^))|n„|^i(Tr)- (42) 

As Sn tends to 0, inequality (|4T]) becomes impossible and proves that Supp(no + uq) has no accumulation points in 
T^. It follows that Uq + no is a finite sum of positive Dirac masses. □ 

Remark 4.2 More precisely, we can get an estimate of the number of corners in each connected component I of 

#{Supp(n[; + no) n /} < ^ + 2 where A'-^ := + 

(C = appears in Lemma |4~T]) . Indeed, let us consider three consecutive Dirac masses 61,62, 63 in /. Then 

• if /3 > 0, 7 > 0, we have 

(ft ft]^-'> -7 + + 4a/? 

• if /3 = 7 = 0, then we have a contradiction, that is to say Uq + no is the some of at most two Dirac masses /. 

To prove this estimate, we define v G Hq{6i,6^) satisfying v" + v = 5q^ in (6'i, ^3), n = in T \ (^i, ^3). In T, 
the measure v" + n is supported in {^i, 62,6^}, and since these points are in Supp(ng + no), and [^i, ^3] C Tj„, 
no + tv is admissible for small The second order optimality condition and then the assumption (fT6l ) together 
with Lemma |4~T] lead to 

< /'(no)(n,n) < -a|w|Hi(T) + 7^1^1(1) Iblk'^m + /^H-fllH^m 
< [-a + C^X + C^|3X^)\v\l^^Jy 

where X = {6^ — 6iY~^ , which implies (1431 ) when j3 is positive, and gives a contradiction if /? = 7 = 0. 

Remark 4.3 When one uses the parametrization of convex sets by the gauge function n, is a polygon if and only 
if u" + n is a sum of Dirac masses. With the support function (see Section [3311, one has the same characterization. 
Therefore, the conclusion is the same if we work with the optimization problem (l39l ). Estimate (1431 ) remains valid. 
However, 9i is no longer the polar angle of a corner of the shape, but is the angle of the normal vectors to the 
successive segments of the polygonal boundary of the shape. □ 

As in Section 12.2.11 one can also handle problem with the equality constraint. 
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Proof of Theorem lllll 

We now need an abstract result for second order optimality conditions. Adapting |[T4l Proposition 3.3] similarly to 
the first order condition given in Proposition 13. II (this explains the assumption /(uq) € (C'^(T))'), we get that there 
exist Co e W^'°°{T) nonnegative, /i G M"^ such that 

{Co = on Supp(tiQ + uq) and 
(44) 
Vug W^'^^iJin), j'{uo)v + n • m'{uo)v = (Co + Co,v){w^^^yxW^'^- 

Furthermore, for all v G H^{Tin) such that 3A G M, witht;" + -y > X{uq+uo), and {Co + Co,v) - iJ.-m'{uo){v) = 0, 

j"{uo){v,v) + • m"{uo){v,v) > 0. (45) 



Then we proceed as in the proof of [14, Theorem 2.1]. Compared to the first step of the proof of Theorem 12.91 we 
add one degree of freedom introducing 4 functions Vn,i on a partition of (0, and we look for Xn,i,i = 1 ... 4 
such that Vn = J2i=i,4 ^n,iVn,i Satisfies 

<(0+) = v'^{e-) = fi ■ m'{uo)vn = 0. 

Such a choice of A„^i is always possible as A„^i satisfy three linear equations. Moreover, Vn is not zero and using 
(|44l) . we get Jj Vn{Co + Co) = 0' which imphes 

= /(no)(u.„) = / VniCo + Co) = M • rn'{uo){vn)- 



As v'^ + Vn > X{uq + uq) for a <^ 0, it follows that Vn is eligible for the second order necessary condition 
Then, it follows 



< j"{uo){Vn,Vn) + fl ■ m" {uo){Vn, Vn) < -a\Vn\H^T) + TI ^^n l/fi (T) ll^n || (T) + (/3 + || ) II 
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As n tends to oo, the inequahty < j"{uo){vn, Vn) + /i • rn"{uo){vn,Vn) becomes impossible and this concludes 
the proof. □ 

Remark 4.4 An estimate similar to the one in Remark 14.21 is not straightforward anymore, since the Lagrange 
multiplier ^ is unknown. 

4.2 Examples 

We analyze the same examples as in Section [l!2l with —P instead of P: 
Example 4.5 (Negative perimeter penalization) One can study 

min{J(0) := F Ef{n), Xi{n)) - P{n) / n convex, Di C C -D2} (46) 

where F : (0, +00) x (—00, 0) x (0, +00) — )• M is C^, / G H'^(M?), and Di, D2 are bounded open sets. We can 
prove the following. 

Proposition 4.6 Ifi^o an optimal set for the problem (I46I ). then each connected component of the free boundary 
BVLq \ {dDi U dD2) is polygonal. 
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Proof. The proof is a direct consequence of Corollary 12.131 and of the estimates given in Section 14.3.21 Indeed, 
Proposition 14. 1 1 I for Ef{Q,), the similar- result for Ai (See Section l4.3.3l) and Proposition 14. lOl for the volume, imply 

|r"(0)| < C\\v\\l,,,^^^^^, 

where r{t) = F{\^t\, Ef{Qt), Xi{^t)), = ^uo+tv and e E (0, |). Next, the estimate for the perimeter in 
Proposition 14. lOl provides the concavity condition. □ 

Remark 4.7 As in Remark 13. 8[ if we consider problems of type (l46l ) where the constraint -Di C C -D2 can be 
dropped, then the solution is a polygon. 

Example 4.8 (Volume constraint and negative perimeter penalization) We can also consider a similar problem 
with a volume constraint: 

mm{J{n) := F{Ef{n), Xi{n)) - P{n) I convex, and = Vq} (47) 

where Vq G (0, +00). Again, Corollary |TT3] applies and leads to the fact that any optimal shape of ( |47] ) is a polygon. 

Example 4.9 [Perimeter constraint] We consider again a problem with a perimeter constraint, as in Example 13.101 

mm{J{n) := F {\n\, Ef{n), Xi{n)) / n convex, and P{n) = Pq} (48) 

where Pq G (0, +00). The optimality conditions are written for £'j(0), Ai(ri)) + iiP{Q,), where ^ is a 

Lagrange multiplier for the constraint P{0.) = Pq, so if we prove that ^ < 0, then the strategy of this section 
applies, and we get that any optimal shape is polygonal. 

4.3 Computations and estimates of second order shape derivatives 
4.3.1 Volume and perimeter 

Let a{u), p{u) be the area and perimeter functional, see (1261 ). 

Proposition 4.10 Let < u G W^'°°{T). Then a and p are twice dijferentiable around u in W^'^{T) and there 
exists some real numbers /3i , /32 , /^s , 7 and a > (depending on u) such that, \/v G W^''^{T) 



(49) 



\a"{u){v,v)\ < I3i\\v\\l2^j^ 

O^Hh^j) -l\v\miT)\\v\\L^T) - /32||^^|li2(Tr) <p"{u){v,v) < /Ss || u|| (-j,) 

Proof. This is done by easy computations, using formulas of Section [3.3.1l □ 
4.3.2 The Dirichlet energy - Torsional rigidity 

We now analyze the second order derivative of e/(n) = Ef{Qu) introduced in Section 1331 The main result is the 
following. 

Proposition 4.11 Assume := Quo> ^0 > 0, Uq + > 0, / G Hf^^{M?). Then ej is in a neighborhood of uq 
(in W^'^{T)). Furthermore, there exist /3i, /32 positive such that, for all v G VF"^'°°(T), 

\e'^{uo)v\ < I3i\\v\\l2(j), (50) 
\e'}{uo)iv,v)\ < /32(lbll^i/2(Tr) + ll^lli-(T))- (51) 
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The differentiability of e/ and the estimate (l50l ) follow easily from Proposition 13.131 The estimate (ISTl) is easy 
to prove when working with smooth sets and one can then even drop the term. However, this result is more 
difficult for a general convex set and the rest of this section is devoted to its proof. 

Let V be given as in Proposition 14. Ill To prove the estimate (|5TI ). it is appropriate to consider a transformation ^ 
such that 

C G C\i-r^, r?), 1^1'°°(]R2, ]R2)), r? G (0, 1), m = f " —] e'' on dUo- (52) 

\Uo + tv UqJ 

Then, we will differentiate twice t G {—ViV) ~^ ^(0 = E{Q,uo+tv)- We will use the following notation and 
identities: 

nt:=nuo+tv, Ut:=Un^^+,,, Ut:=Uto{I + i{t)), e{t) := E{nuo+tv). (53) 
Note that e{t) = ef{uo + tv) = <f/,f7o ^^id we have 

e'{0)=e'f{uo){v) = 4f,„(0)(e'(0)), (54) 
e"{0) = e'}{uo){v,v) = S'luoim' {0), ^ (0)) + £'f,no im"m. (55) 

In the smooth case, e"(0) can be written in terms of boundary integrals, which involve in particular the boundary 
trace of D'^Uq and VC/q. These terms are not well defined in the non-smooth setting (even in the case JIq convex). 
To overcome this difficulty, our strategy will be to write all non-smooth terms of e"(0) as "interior" integrals in Qq. 

Estimate of e"(0): Note that we have proven in Section [33] that e/ is if / G Hi^^{R'^)) (so, e is C^). We remind 
the following classical formulation of e"(0) 

Lemma 4.12 Let f G Hl^{M.'^) and ^ G C^{R; VF^'°°(M2; M^)) near 0. Then we have 

e"(o) = -\{l m + 1 Miem ■ >^o)) , (56) 

where Uq := U^{uo){v, v) and Uq is defined by 

C/^' := ^7^' - (2Vf/^ e'(0) + C'(0) • Z?'f/o • e'(0) + Vf/o • e"(0)) in Qq, (57) 

and satisfies 

Uq G L^i^o), AC/q = in P'(Oo). (58) 
Proof. Differentiating (|37] ) at u = uq (see Corollaire 5.2.5, ifTOl ) and then using (l34l ) gives 

e"(0) = ({>^'-2V^^C'(0)+e'(0)-I)'C/o-C'(0)+V^o-(2VC'(0)-e'(0)-C"(0)))/, 

~l [ UU{^'{0)-uo). (59) 

After replacing = - VC/q • C'(0), <l59ll gives (l56ll. 

Clearly Uq G L^(r2o)- To prove that A[/q = we differentiate (1351 ) atu = uq and use (l34l ). Then we obtain 

^ ({7^' - 2VC/^ • C'(0) + e'(0) • D^Uo ■ e'(0) + VUo ■ (2VC'(0) • ^'(0) - ^"(0))) Aip = 0. 
Replacing Uq as given by (|3T]) gives J^^ UqAip = 0, which proves (l58l) . □ 
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Proof of Proposition 14.1 1[ 

We will often write , ^" for C(0), C'(O), ^"(0). Let us rewrite §B in the form e"(0) = ^{h + h)- The second 
term I2 is easy to estimate: from (|3TI) we have 

h := - I m{e-i^o)=[ fduoUoie-i^o?<CU'\\l2^9n), C = C{\\fU^^n,),\\Uo\\w^,o^inopO) 

The first term Ii = J^^ UqAUq requires more investigation. To go around the non regularity of r2o> we introduce 

Uo = Ui-U2,U^€H^{no),-AUi=f+, -AU2 = r, C/^ > on ^^q. 

Recall that Ui G W^'^^iilo) n H'^i^o). We will compute on the level sets Ql := {x G Qq, Ui{x) > e} (only on 
one of them if = or /" = 0). Indeed, by Sard's theorem, the 0* are at least for a.e. e. By strict positivity 
of Ui, lim^^o = 1^0' so that 

h = lim [ UqAUi - [ U0AU2. 
^-^oJni Jni 

Note that C/^, G C^^i^o) and as / G HI^{R'^) we have Uq G Hf^^{QQ). We obtain 

=: II + 11 + 11 + 11 (61) 

For the term If, we have 

If = I U'^AUi + VU'^-VUi^ I U'^AUi + VU'^-VUi= [ U'^d,,Ui = ^. (62) 

To deal with If and If, we will need the following generalized formula of integration by parts. 

Lemma 4.13 Let Q be a open set, U G Wl''^{n)r\H^{il), V G H^{n)r\{AV G L'^{^q}, g G M^). 
Then 

J := [ dyU{g-W)= [ V{VU ■ g) -VV + {VU ■ g)AV -V{V^U ■ g) -V^V, (63) 
Jan Jn 

where the operator _L acts on a vector and is defined by """(ai, 02) = (— ^i)- a consequence 

\J\ < ||VC/||i«=(n)||sr||i2(Q)||Ay||i2(Q) +2{||y||^i(Q) [||C/||//2{n)||s'llL-{Q) + l|V[/||L^(n)||V5||i2(f^)] .} (64) 

Proof. If V is the exterior normal unit vector to dO. and r =^v the unit tangent vector, then, for Lp G H^{^) and 
a = (oi, 02), using that VC/ • r = on dQ.Q, we have 

^r = -z/, {a-v)dyU = a-VU, {a ■ T)dyU = a -V^U, t -Vip = -v -V^^. 

Then we obtain 

/ d,U{g-VV) = [ {VU ■ i^){g ■ i^){VV ■ u) + {VU ■ i^){g ■ t){VV ■ t) 
Jdci Jan 

= / ( (VC/ • g)W — {V'^'U ■ g)'V^V] ■ v (apply divergence theorem to both terms) 

ho. ^ ^ 

JO. 
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which proves (l63l) because V • V"*- = 0. The estimate (l64l ) follows. 
End of the proof of Lemma 14.131 



□ 



We apply Lemma l4.13l on Q = Ql to estimate If, If in (|6TI ). For /I, we choose U = Ui — e,V = U^^g = ^'(0) 
(recall that /^U'q = 0) and for /|, we choose U = Ui - €,V = Vj = djUi,g = gj = Cj(0)C'(0),i = 1,2: 
here — AV, = djf'^ or djf~. Next, we apply the estimate (l64l) to each of these choices and we are obviously led 
to estimates independent of e. For J|, we make a direct easy estimate. Together also with (l60l) and using Young 
inequality we obtain: 

|e"(0)| < C(||VC/^||i2(f,„) + ||Ve'||i2(f,„) + U'Wl^^no) + U'WUdno) + WUHano)), (65) 

where C = C{\\f\\Lo.^n,,)nm{no)^ \\Ui\\w^,^(no)^ \\Ui\\H^no)^ ^ = 0' l'^)- (66) 

Now, let us write the estimate (1651 ) in terms of v. First note that if a, /3 G i/^/^(9r2o) H L°°{di}) then a/? € 
Fi/2(517o)nL°°(5fio)and 

ll"/5|lHi/2(ano)nL°°(9f7o) - ^ll«llHi/2(af7o)nL°°(ano)ll/^llHi/2(ano)nL°°(9t7o)- 

(using the easy fact that H^{Uo)f^L°^{0,o) is an algebra, and that the H^^^ {dO,o)-norm is equivalent to the H^{Qo)- 
norm of the harmonic extension in Jlo)- Also, we point out that the transformation tjj = ilj{r,9) := ^J^g^ e^^ is 

bi-Lipschitz near- T and ^(T) = 90o- Then 7 G H^/^{dn) if and only if 7 o g H^/'^{T), and their if V2 .norms 
are equivalent. 

Let us remind that, according to the choice of ^ in (l52l) . we have ^'(0) = — -^e*^, ^"(0) = 2^ on d^Q. Then we 
obtain, with the same dependence of the various constants C as in ( 1661 ) 

l|V^ollL2(f7o) - ^1^' ■ '^^01^1/2(0!^) - '^ll^'llHi/2(aQ)nL°°{no) - '^ll^llHi/2(T)nL°°(T)' (6"^) 

ll?'lli,°°{eno) < C'lbllL°°(T), ll'?'lli2(eno) + ll^"llii(9f^o) ^ C'll^lli2{T)- (68) 
All these estimates are valid for all choices of ^ as in ( |52l ). Let 

Given w G W, let us choose ^(t) := + t{w — C'(0)), where ( is the l^^'°°-extension as given in ( [28] ). namely 

Vuo(6') / \uo{9) + tv{e) uo{9)J 

with 7/ = (resp. t] = 1) in a neighborhood of the origin (resp. of d^o)- Then, ^ is as in (l52l ) and ^'(0) = ij;. 
Therefore, the estimate (l65l) together with (l67l) . (1681 ) leads to 

Ww G W, |e"(0)| < C [\\Vw\\l,^^^) + lbll^i/2(Tr)nL-(T) + ll«^lli-(f7o)) • (69) 

Let us introduce 

wo£H\no), Au;o = OonOo, (u^o)|ano = "^e^^ [or ^i;o - C'(0) G i^o^^^o)] • 
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Let now 5„ be a sequence of C^(r2o)-functions converging to wq — C'(0) in Hq{Uq) and let Wn '■= inf{5„ + 
C'(0), ||ifo||L°°(no)}- Then, Wn ^ W and converges in if^(ilo) to wq. Applying (|69l ) with Wn in place of w and 
passing to the limit yields: 

|e"(0)| < C [WVwoWh^n,) + Mm/^T)nL^m + lho|lioo(f,„)) . (70) 
But, since wq is harmonic, 

o) — ll''^0||L°°(9r2o) ^IIl°°(t)- 

Finally, the estimate (TTOI ) leads to 



2"(0)| < C7||7;||^i/2(Tp)n^^(Tr). 



□ 



4.3.3 First eigenvalue of the Laplace operator with Dirichlet boundary conditions 

The estimate of Proposition 14. 11 1 also holds for Ai(r2M), the first Laplace eigenvalue (see Section [331 ). namely 

|/i(0)| < C\\vfH^^\T) n L°°(T), (71) 
where li{t) = Xi{Quo+tv)- As the computations are very similar, we will only sketch the proof. 

Proof of (ItB. 

As for e/, for v € W^'°°{T) fixed and |t| small we consider li{t) := Xi{^t) and Ut, the first eigenvalue and the 
corresponding eigenf unction of — A in := ^luo+tv As in Lemma 14.121 we can show that 



^i(O) = - [ U'^^U'^ + UoAU;; =: h + h. (72) 

Here U'q and Uq satisfy 

l^ol' + Wo =0, 



where li = /i(0), /'^ = /'^(O), I'l = /'/(O). Then considering Q.^ = {x ^ U,q,Uq > e} as in the proof of Proposition 
14. llK note that J/q > on Q.q here), we have : 

h = f KihUl, + l[Uo) = h [ \U(,\\ (73) 



-Au;, = hu;, + i[Uo 


in Oq, 


Uo = -e'(0) • VUo on 




/ 






= hU'^ + 2l[U'Q + l'iUQ in 




/l 


U'^ - (2C'(0) • V% + e'(o) 


• D^U^ ■ 


e'(0)+e"(0)-Vf/o) in 







2 



- lim /" UqAUq = - lim /" eAC/o + (Uq - e)AUQ = - lim /" {Uq - e)AUo 
lim / Udi-AUo) - lim / ([/q - e)a,.C/^' - C/^'9.,f/o 

-h [ K\^ + lirn [ U'^d^^Uo. (74) 
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Combining (1721 ) with the last two equalities gives 

r/(o) = lim / u;^dM- 

Jan, 

Then we proceed exactly as in Proposition 14. 11 [ and obtain for l'{{0) an estimate exactly similar to (|65] ). 

Next, we prove that WVU^^h^no) < <^ll^'(0)||//i(no)- As = (j^ - ^'(0) • VC/q it is enough to prove 
l|VC/ollL2(no) < C'IIC'(0)||//i(no)- One can verify that % satisfies 

U'q £ H^{no) n H\no), Af/^ + hU;^ = 2trace([V^'] • [D^Uq]) - l[Uo, [ U'^Uq = 0. 

Using the convexity of VLq and Fredholm alternative theorem, we can prove that the operator 

y G (/7o^(!^o)ni72(f]Q))\span{C/o}^ Ay + /iT/GL2(!^o)n{/i, f /if/o = 0}, 

defines an isomorphism (see for example lEl), which together with the formula for /'i(O) provides the required 
estimate for Uq. Therefore, as for e"(0), for all ^ as in (l52l ). we have 

//i/2(T)nL°°(T) °°(f^o) ) • 

Then we complete the proof as in Proposition 14.111 □ 

5 Remarks and perspectives 
5.1 Localization of our two approaches 

As explained in the introduction, the approaches leading to our two families of results are very "local" with respect 
to the boundary of the optimal shape. Indeed, each proof uses test functions v G whose support may be as 

small as we want and only covers the portion of the boundary that we want to analyze. To show how this can be 
exploited, we give now -without proof- , an example of a result which can be reached by the same two methods 
when applied locally. 

Let us consider the following optimization problem where G : (0, g') G T x M x M — )• M is assumed to be of 
class and o, 6 G (0, oo): 

' uq G 1^^'°°(T), j{uQ) = min{j(n), n" + n > 0, a<u<b}, 

r (75) 

where j(u)= / G {0 , u{e) , u' {0)) dO . 
Jt 

We define as in ([T3l) and we introduce the partition Tf„ = T+ U Tq U T_ where 

T+ := {e G Tin ; G^(^) G (0,oo)}, (recall G^) = Ggg {0,uo{0)Mo{0+) + (1 " 0^o(^~)) dt, 
T_ := {6 G Ti„, [G,g {0, uo{0), u'o{0~)) , G^, {6, uo{e),u'o{e+))] C (-oo, 0)}, 
To :=Ti„\(T+UT„). 
Then 

(i) T+ is open and u'^ G L^^{T+), so that uo G W^^^{T+), 

(ii) There is no accumulation point of Supp(uQ + uq) in the open set T_; in other words, [9 G T_ — dQ,uo (^)] is 
locally polygonal. 

The situation on Tq requires a complementary study specific to each functional. 
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5.2 Very singular optimal shapes 

In this paper, we gave some sufficient conditions on the shape functional so that an optimal shape be smooth 
or polygonal. But there exist convex sets which are not of this type, and in a certain sense have "intermediate 
regularity". Namely, there are convex sets which are singular in the sense that they do not have corners (they are 
C^), but their curvature is zero almost everywhere. As an example, one may consider any convex set such that 
u" + n is a Radon measure, without mass, but singular with respect the Lebesgue measure. 

Let us mention a shape optimization problem whose solution is neither regular nor polygonal (see IITtI for an 
analysis of this problem). Let be a convex set, Vq = \^o\, Pq = P{^o) and D = (Oo)t = {x E M?, d{x, Qq) < 
T}. Then Theorem 8 in Inl states that: 

J{no) = min{ J(0) / nc D convex such that = Pq, \n\ = Vq}, (76) 

where J is the distance functional: 

J{n) := [ d{x,n)dx. 
Jd 

Since Q is any convex set, one cannot expect any geometrical property for a minimizer of (1761 ) without extra condi- 
tions on D, Vq and Pq. Remark also that the box D = (r2o)r is C^'^ here. 

5.3 Problem without perimeter 

An interesting problem, which has not been analyzed in this paper, is the following (we use the notation of Section 
IH: 

max{£'/(r^), \Q\ = Vq, O convex C D}. (77) 

It is easy to prove the existence of an optimal shape r^o- In this situation, we expect the term Ef{Q) to be leading 
over \ fl\ (whereas the perimeter was the stronger term in the examples solved in this paper). So we are naturally led 
to the following question : do there exist a > 0, /3, 7 > such that 

yveW^'^^iT), e"{0)>a\v\l,^,-^\v\Hi/2\\v\\L2-P\\v\\h'^ (78) 

A consequence of such an estimate, would be that any solution of dTT] ) is locally polygonal inside D (the same 
strategy as in the proof of Theorem |2]9] would provide the result, we just need to adapt Lemma |4TT] to ff^/^-norms). 
It is easy to prove that (1781 ) holds if v is supported by a subset of T which parametrizes a strictly convex part 
of dQo. Therefore, with the same proof as for Theorem 12.91 we are in position to deduce that dQo H D is nowhere 
with a positive curvature. But it is not cleai^ whether estimate (1781 ) remains valid in a more general situation and, 
consequently, whether OUq n is a polygon or not. 
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